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A passive neutron multiplicity measurement device, FH-NCM/S1, based on field-programmable gate arrays 
(FPGAs), is developed specifically for measuring the mass of plutonium-240 (**°Pu) in mixed oxide fuel. FH- 
NCM/S1 adopts an integrated approach, combining the shift-register analysis mode with the pulse-position 
timestamp mode using an FPGA. The optimal effective length of the *He neutron detector was determined to 
be 30 cm, and the thickness of the graphite reflector was ascertained to be 15 cm through MCNP simulations. 
After fabricating the device, calibration measurements were performed using a *’Cf neutron source; a detec- 
tion efficiency of 43.07 % and detector die-away time of 55.79 us were observed. Nine samples of plutonium 
oxide were measured under identical conditions using the FH-NCM/S1 in shift-register analysis mode and a 
plutonium-waste multiplicity counter. The obtained double rates underwent corrections for detection efficiency 
(£) and double gate fraction (fa), resulting in corrected double rates (D.), which were used to validate the ac- 
curacy of the shift-register analysis mode. Furthermore, the device exhibited fluctuations in the measurement 
results, and within a single 20-s measurement, these fluctuations remained below 10%. After 30 cycles, the 
relative error in the mass of "Pu was less than 5 %. Finally, correlation calculations confirmed the robust 
consistency of both measurement modes. This study holds specific significance for the subsequent design and 


development of neutron multiplicity devices. 
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I. INTRODUCTION 


Mixed oxide (MOX) fuel constitutes approximately 5 % of 
the global inventory of light-water reactor fuel, with an annual 
utilization rate of approximately 2000 tons. Forecasts suggest 
that MOX fuel is likely to play a pivotal role in breeder reac- 
tor initiatives globally. The uranium/plutonium equilibrium 
within MOX fuel is a pivotal parameter that requires meticu- 
lous scrutiny. Both the uranium and plutonium compositions 
within the fuel and segregated uranium and plutonium con- 
stituents must be determined to achieve a comprehensive and 
closed uranium/plutonium mass equilibrium [24—26]. 

Nondestructive assay (NDA) methodologies and appara- 
tuses for determining uranium and plutonium contents rep- 
resent a focal point in nuclear-fuel assessment. This in- 
volves the utilization of techniques, such as detecting ura- 
nium and plutonium isotope-decay processes, particularly 
through the observation of y radiation emanating from nu- 
clear fuel [5, 6, 9]. For example, Canberra’s Tomographic 
Gamma Scanner (TGS) is a representative product [8] capa- 
ble of providing comprehensive quantitative gamma analy- 
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sis for fissile materials or transuranic radionuclides. How- 
ever, its applicability is constrained by factors such as pro- 
longed measurement duration, intricacies in equipment op- 
eration, associated costs, specialized skill requirements, and 
the exclusive focus on evaluating plutonium quality in minute 
uranium-plutonium amalgams or plutonium oxides, render- 
ing this methodology prohibitively expensive. 

Quantifying plutonium content based solely on y-ray [27] 
intensity is compromised by the inherent self-absorption ten- 
dencies of the sample. Notably, uranium and plutonium iso- 
topes not only emit y rays, but also exhibit pronounced mul- 
tiplicity [1-3] characteristics through the release of robust 
spontaneous fission neutrons. Capitalizing on this principle, 
the High-Level Neutron Coincidence Counter was developed 
in 1979 by the Los Alamos National Laboratory (LANL) in 
the United States to gauge substantial plutonium samples. 
The methodology was adopted by the International Atomic 
Energy Agency for global plutonium material accounting [4]. 

Canberra innovatively introduced a suite of products based 
on the spontaneous fission of “Pu and the concurrent emis- 
sion of multiple fast neutrons. These products leverage three- 
ring, four-ring, and five-ring multiple neutron detectors for 
coincident measurements, achieving a detection limit of less 
than 10 nCi/g with a sample volume of 208 L. 

In 2009, the Chinese Academy of Atomic Energy Sciences 
pioneered the development of a neutron-measurement appa- 
ratus for nuclear waste encapsulated in drums. This device 
quantified the mass of Pu by assessing the neutrons gener- 
ated through spontaneous fission. Remarkably, the detection 
efficiency of the device was 30.6 %, with a lower limit of de- 
tection below 100 mg of the effective mass of 74°Pu. 

Owing to the scarcity of *He gas, the focal point of in- 
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Fig. 1. Diagram of FH-NCM/S1 passive neutron multiplicity device. 


strumental research and development (R&D) underwent a 
paradigm shift towards innovation in novel neutron detectors 
post-2010. Multiple R&D teams in cognate domains have pi- 
oneered the development of a suite of coincidence measure- 
ment and multiplicity devices. These innovations are based 
on boron-coated proportional counters [17], liquid scintilla- 
tor detectors, and LiF/ZnS-based scintillator neutron coinci- 
dence counters [20]. However, these advancements are ac- 
companied by challenges, such as suboptimal detection effi- 
ciency [5] and protracted neutron decay times, necessitating 
extended measurement durations to ensure the precision of 
each assessment [20, 22, 23]. Therefore, despite these en- 
deavors, *He neutron detectors are preferred for applications 
within the industrial sphere. 


Two methodologies are currently dominant in the field of 
data processing across various devices. The first is based on 
shift registers, as exemplified by the JSR-15 Handheld Multi- 
plicity Register crafted by the LANL. However, data acquisi- 
tion entails the amalgamation of data outputs in list mode, a 
paradigm embodied by the Advanced List Mode [15] Multi- 
plicity Module [19], also manufactured by the LANL. 


This study investigates the requirements for a nondestruc- 
tive assay (NDA) of 74°Pu quality in MOX fuel, presenting 
the conceptualization, assembly, and evaluation of an inno- 
vative plutonium quality-assessment apparatus, FH-NCM/S1. 
To optimize the instrument cost while adhering to the user- 
specified technical parameters, a Monte Carlo simulation is 
employed to determine the effective length of the primary 
component, 3He, and dimensions of the graphite reflector 
layer. 


In pursuit of diversity in analysis circuits, a holistic strat- 
egy is adopted by integrating the multiplicity register with 
the list-mode multiplicity method using field-programmable 
gate arrays (FPGAs). This approach offers users an expan- 
sive array of data-acquisition methods. Each *He neutron de- 
tector is accompanied by a preamplifier circuit, wherein the 
nuclear signal generated undergoes amplification and shap- 
ing and is transformed into a square-wave signal with a pulse 
width ranging from 2 to 3 microseconds. After scanning the 


pulse signal via the FPGA, the multichannel outcomes are 
consolidated into a singular digital signal after noise reduc- 
tion. 

The FH-NCM/S1 is designed for **°Pu-quality measure- 
ment and underwent rigorous analysis through experimental 
assessments using a 7>*Cf neutron source and standard sam- 
ples. The rationality and precision of the instrument are au- 
thenticated by cross-referencing measurements of identical 
standard samples using a traditional plutonium scrap multi- 
plicity counter (PSMC). 


Hf. DEVICE DESIGN OPTIMIZATION AND RELATED 
PARAMETER INDICATORS 


A. Overall Design of the Device 


The passive neutron multiplicity measurement device con- 
sists of a relatively fixed structure, including a sample cham- 
ber, top- and bottom-end plugs, neutron detectors, modera- 
tors, and electronic circuits. Fig. 1 shows a schematic of the 
passive neutron multiplicity measurement device designed to 
measure the mass of 74°Pu in special materials, named FH- 
NCM/S1. 

The left-hand side of Fig. 1 outlines the primary structure 
of the device. The sample chamber is positioned at the cen- 
ter and adheres to the design specifications, which assume 
a cylindrical configuration with a bottom diameter of 20 cm 
and height of 20 cm. Two end plugs primarily composed of 
graphite were placed at the top and bottom of the sample cav- 
ity. These end plugs served as reflective layers that decreased 
the escape of fast neutrons generated by sample fission from 
the axial direction of the sample cavity. The *He neutron- 
detector array system is shown in blue. This array was ar- 
ranged in high-density polyethylene in four concentric rings 
according to a predetermined pattern to detect spontaneous 
fission neutrons within the sample. 

The overall detection efficiency of the device is a crucial 
indicator of its design process. Both the effective length of 
>He in the detector and the thickness of the graphite reflector 
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Fig. 2. FPGA signal acquisition and analysis process (shift-register analysis mode). 


layer on the top plug play pivotal roles in determining this ef- 
ficiency. Monte Carlo software simulations are subsequently 
used to establish these parameters, ensuring that the device 
design meets the criterion of detection efficiency exceeding 
40 %, while also reflecting its economic feasibility. 

To prevent the diffusion of thermal neutrons back into the 
sample chamber, a 1-mm-thick cadmium layer was incorpo- 
rated on the exterior of the sample chamber to absorb thermal 
neutrons, as illustrated in Fig. 1 (1) and (2). Additionally, 
1-mm-thick cadmium layers were implemented on the outer 
side of the inner slowing-down body and on the inner side of 
the second slowing-down body. The cadmium layer on the 
outer side of the inner slowing-down body absorbs thermal 
neutrons leaking from the sample that are undetected by the 
detector, as depicted in Fig. 1 (3). Moreover, the cadmium 
layer on the inner side of the second slowing-down body ab- 
sorbs thermal neutrons entering the device from the outside, 
thereby reducing the impact of the environmental neutron 
background on the device performance, as shown in Fig. 1 
(4). 

The signal acquisition, analysis, and processing module of 
the FH-NCM/S1 is executed through an FPGA-based embed- 
ded system capable of enabling high-speed digital signal pro- 
cessing and communication functionalities. A 50-MHz clock 
crystal oscillator ensured a measurement-time accuracy of 20 
ns. The device incorporates both the shift-register analysis 
mode (signal analysis and processing are completed through 
the shift registers and up-down counters), as illustrated by the 
green block in Fig. 1, and the list mode (referred to as the 
pulse-position timestamp mode in the text), depicted by the 
blue block in Fig. 1. Users have the flexibility to choose a 
suitable measurement method through the software interface, 
enabling the device to adapt to diverse measurement scenar- 
10S. 

The traditional shift-register analysis mode consolidates 
multiple signals into a unified form, channeling the digital 
signal generated by the pulse time series into a shift register. 


The distributions of R+A (real count and accidental count) 
and A (accidental count), which are crucial for multiplicity 
analysis, are derived by sequentially inputting the digital sig- 
nal into distinct shift registers, such as the pre-delay time, 
coincidence gate width (or gate width), and long delay time, 
while regulating the output. For ease of portability, the sig- 
nals from the eight *He neutron detectors are typically amal- 
gamated into a singular signal via an OR gate [4] and subse- 
quently routed to JSR-15 for analysis. 

However, in the design of FH-NCM/S1, careful consider- 
ation is given to the self-check predicament inherent in the 
detectors. The device employs an FPGA to directly scan 80 
signals from the *He neutron detector, enabling multichannel 
signal acquisition and simultaneous recording of the origin of 
the signal (ID of the detector). The analysis method is shown 
in the green block in Fig. 1. 

The FPGA now encompasses the pulse-position timestamp 
mode, enabling not only direct industrial multiplicity anal- 
ysis by the shift-register analysis mode, but also providing 
intricate pulse-time information. The output mode of Atp, 
representing the time interval between two recorded pulses in 
the pulse-position timestamp mode, is more suitable for jitter 
filtering and noise removal. Simultaneously, this mode can 
retain the original signal-measurement content for subsequent 
viewing and offline analysis. The output result is illustrated 
in step 2 of Fig. 1. 


B. Principle of FPGA Signal Acquisition and Analysis 


In the shift-register analysis mode, the front end of the 
data-acquisition and analysis device comprises several OR 
gate circuits. These circuits primarily eliminate randomiza- 
tion from multiple signals, consolidate them into a unified 
output, and facilitate subsequent connections with the acqui- 
sition and analysis device. However, this processing method 
requires a high level of consistency in the output-pulse signal 


from each *He neutron detector and does not log the origin 
of the pulse signal (i.e., which *He neutron detector generates 
the signal). This poses a challenge in identifying and address- 
ing the issues related to the failure of *He neutron detectors. 

To overcome this limitation, the signal-acquisition and 
transmission module of FH-NCM/S1 adopts a direct connec- 
tion method. The output signal of the *He neutron detector 
is directly linked to the FPGA pins, enabling the scanning 
and collection of all the pulses generated by the detectors. In 
this process, the analog pulse signal from the *He neutron de- 
tector undergoes preamplification, shaping, and amplification 
with a determined amplitude threshold. The resulting dig- 
ital pulse signal is then directly fed into a digital-processing 
chip FPGA, where the pulse is recorded. The FPGA precisely 
records the position (Detector ID) and time of the acquired 
pulse front, ensuring a minimum resolution of 20 ns. 

Regarding the analysis principle, FPGA functions are 
based on the acquisition method of traditional multiple-shift 
registers. The pulse time-stamp captures the collected signal 
as a time series, subsequently inputting it sequentially into 
the corresponding "coincidence" gate for shifting (as depicted 
in step 1). Within both the (R+A) and A gates, an up-down 
counter is programmed to record the data in real time in the 
gate (as illustrated in step 2). Upon completion of the data- 
shifting process, the results recorded by the reversible coun- 
ters in the (R+A) gate and the A gate are the output (as indi- 
cated in step 3). A visual representation of the FPGA setup is 
shown in Fig. 2. 

The FPGA accumulates statistical information from the 
(R+A) and A gates, which is then transmitted to the host com- 
puter. The host-computer software calculates the single/dou- 
ble/triple rate (S/D/T) using Eq. 1, Eq. 2, and Eq. 3 in subse- 
quent analysis. The distribution in the (R+A) gate is denoted 
as P(n), the distribution in the A gate is denoted as Q(n), and 
‘t represents the measurement time [4]. 
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Furthermore, based on Eq. | and Eq. 2, the single rate rep- 
resents the total neutron counting rate detected by the detector 
array, whereas the double rate represents the real count rate, 
that is, the coincident count rate. Compared to S, D more ac- 
curately reflects the rationality of the algorithms used in the 
device and the variations in various setup parameters. There- 
fore, it is often employed as a benchmark to analyze the ra- 
tionality of a device. 


In addition to the direct output of (R+A) and distribution of 
A, the shift-register analysis mode can be switched to the po- 
sition timestamp mode using the host computer software. The 
At value between adjacent pulse signals can be directly out- 
put and subsequently integrated into a comprehensive pulse 
time series using the upper computer software for analysis. 
This approach facilitates the acquisition of a complete mea- 
surement response for the detector. It expedites the identifi- 
cation of any detector issues and enables prompt adjustments 
such as detector replacement or the fine-tuning of denoising 
algorithms. Specific details are discussed in Section IV. B. 


C. Device Parameters and Evaluation Methods 
I. Solution of Neutron Multiplicity Counting Method 


The neutron multiplicity counting (NMC) method deter- 
mines the mass of *“°Pu in a sample by formulating and 
solving measurement equations. These equations rely on 
the assumption that the sample can be approximated as a 
point-source model situated at the center of the sample cav- 
ity. The measurement equations are derived by solving the 
multiplicity distribution of the source-event neutrons using a 
probability-generating function, as shown in Eq. 4, Eq. 5, and 
Eq. 6. The single rate (S), double rate (D), and triple rate (T) 
for the NMC method are obtained using the following expres- 
sions [2, 4, 29]: 


S = FeMvs, (1 +a) (4) 
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In this equation, in addition to S/D/T, F denotes the sponta- 
neous fission rate of the sample, which represents the overall 
fission level of the sample in units of st. e denotes the de- 
tection efficiency of the instrument for neutrons and M is the 
neutron-leakage multiplication of the sample. The terms Vsn 
and vi, denote the first, second, and third reduced moments of 
spontaneous and induced fission neutron distributions in the 
sample, respectively. œ represents the ratio of the number of 
(a,n) neutrons produced by a reaction to the number of neu- 
trons produced by fission. fy and fı are the double and triple 
gate fractions, respectively. 

Eq. 4, Eq. 5, and Eq. 6 treat the fission-material problem 
as a parameter-estimation challenge and utilize the moment- 
estimation method to determine the neutron-leakage multi- 
plication M. Eq. 7 establishes the relationship between M and 


the constant coefficients a, b, and c in the equation, which can 
be expressed as Eq. 8, Eq. 9, and Eq. 10, respectively [4]. 
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By substituting the neutron-leakage multiplication M of the 
sample into Eq. 11, the spontaneous fission rate of the sample 
can be determined. 
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For the plutonium samples, where the fission rate per unit 
mass of *4°Pu is 473.5 gtst}, the equivalent mass of 74°Pu in 
the sample can be expressed using Eq. 12: 
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If the isotopic abundance of the plutonium is known, the 
calculated plutonium mass, denoted by m., can be obtained 
(Eq. 13), where fn represents the proportion of the different 
isotopes. 
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After calculating the mass (m<) of *°Pu in the sample, the 
relative error of the device measurement is determined by 
comparing it with the actual mass (ma) of 7*°Pu in the ex- 
ample using Eq. 14. 


A= 


Ze Te x 100% (14) 


Ma 


However, before employing the NMC method to ascertain 
the mass of **°Pu in the sample, the parameters such as €, fa, 
and f, must be calibrated. As fa and f, are functions of the 
detector die-away time (T), both £ and 7. 


2. Detection Efficiency 


Because Pu is classified as a military-controlled substance, 
device calibration often utilizes a **Cf neutron source with 
energy levels similar to those of **°Pu fission neutrons. Eq. 15 
presents the calculation formula for the detection efficiency. 
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Among these parameters, Ngerecteq represents the total count 
rate of the detection array comprising 80 *He neutron detec- 
tors and Nemitea represents the neutron yield of the known neu- 
tron source. In practical calibration, a 252Cf source is com- 
monly used to expedite the calibration process. To mitigate 
systematic errors arising from environmental and electronic 
circuit influences, the background impact should be consid- 
ered when calculating Ngetectea. The mean results of multiple 
measurements are used to reduce the impact of accidental er- 
rors. Subsequently, Ngerecteq can be expressed using Eq. 16. 


"detected = No Eg No (16) 


Here, No represents the average active count rate obtained 
from repeated measurements of the instrument when the >*Cf 
source is at the center of the sample chamber, and N, repre- 
sents the average background (BG) count rate. 


3. Detector Die-away Time 


The detector die-away time (7) is another crucial device 
parameter, along with the detection efficiency (£). This signi- 
fies the rate at which the neutrons decay from the same fission 
event in the device. 

To ascertain the detector die-away time (7), one approach 
involves measuring the double rate (D) of a single neutron 
source at various ‘gate widths’ and employing the relationship 
in Eq. 17. The detector die-away time (7) can be determined 
through curve fitting, with Do representing the double rate for 
an infinite gate width. 


D(GW) = Do x f — el- F] (17) 


The curve-fitting method helps alleviate the computational 
challenges stemming from the limited selection of gate-width 
combinations. However, to enhance the precision of the final 
result, more data points must be measured, and the error in 
each measurement result contributes to the overall error. Con- 
sequently, when employing the curve-fitting method, repeated 
measurements must be conducted under the same gate-width 
value to minimize the impact of accidental errors. 

After constructing the device, considering a practical sce- 
nario, each measurement lasts for 20 s, and the process is 
repeated 30 times. The measurement results are recorded and 
the 7 values of the device are calibrated according to Eq. 17. 


Eq. 18 and Eq. 19 represent the double- and triple-gate frac- 
tions, fq and f;, respectively. PD represents the pre-delay time 
in microseconds, and G denotes the gate width in microsec- 
onds. The actual values can also be calibrated by using mul- 
tiple Cf sources with different neutron yields. 
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4. Repeatability of the Measurement Results 


In this study, the standard deviation of a single measure- 
ment is used to represent the measurement repeatability RSD 
of the instrument. The calculation formula is shown in Eq. 20, 
where mi is the calculated mass of the same sample for differ- 
ent measurements, mM is the average value of the sample mass, 
and n is the number of measurements. 


ae (mi—m)? 
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RSD = x 100% (20) 


240Pu samples with various masses were subject to repeated 
measurements. To align with practical conditions, the indi- 
vidual measurement time was established at 20 s, and the 
measurement cycle was repeated 30 times. Following these 
measurements, the RSD values of the results were calculated 
using Eq. 20, which facilitates the analysis of the repeatability 
of the device-measurement outcomes. 


HI. DESIGN PARAMETER DETERMINATION AND 
DEVICE CONSTRUCTION 


A. Simulation and parameter determination 


Based on the design plan, the “He neutron detector-array 
system comprises 80 detectors arranged in a circular pattern 
around the central axis of the sample chamber. The detec- 
tors are organized into four rings, with each ring containing 
19, 25, 18, and 18 detectors. Furthermore, as determined in 
Section II.A, the effective length of 3He and thickness of the 
graphite reflective layer atop the plug are indispensable fac- 
tors that have a profound impact on the capacity to discern 
efficiency. To achieve a detection efficiency above 40% at 
the lowest cost, Monte Carlo software is used to simulate and 
analyze various factors that affect the detection efficiency. 
This meticulous procedure is undertaken with the sole aim 
of deducing the optimal values of the effective length of He 
and the thickness of the graphite reflective layer. The cho- 
sen simulation software is MCNP, with various nuclear reac- 
tion cross-sections referenced from the NDF/B-VUI.0 nuclear 
database [30, 31]. 
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Fig. 3. Detection efficiency of the device varying with the effective 
length of the detector. ((a) Influence of effective length of 3 He neu- 
tron detectors; (b) influence of graphite feflective layer) 


During the simulation, a Cf point source is strategically 
placed at the focal point within the confines of the sample 
chamber. The sample chamber is maintained at a height of 20 
cm, and the effective length fluctuates between 25 and 45 cm. 
Fig. 3 (a) illustrates the relationship between the detection 
efficiency of the device and effective length of the detector. 
From Fig. 3 (a), the detection efficiency increases with an 
appropriate sensor size. When the detector’ size is 30 cm, the 
detection efficiency surpasses 40 %. Therefore, a *He neutron 
detector with an effective length of 30 cm proves to be the 
most cost-effective and efficient choice for a detection array. 


Considering the distribution of the *He neutron detectors, 
the thickness of the graphite reflective layer is progressively 
increased from 0 to 15 cm. Fig. 3 (b) presents the overall 
detection efficiency of the device and trend of the detection 
efficiency per ring as the thickness of the graphite reflective 
layer increases. The results indicate that the detection effi- 
ciency increases with the thickness of the graphite reflector. 
Owing to the space limitations in the sample cavity, the thick- 
ness of the graphite reflector in this study is set to 15 cm. 

Each He neutron detector has a diameter of 2.54 cm and 
an effective gas length of 30 cm. The upper segment incor- 
porates electronic circuits extending approximately 15 cm in 
length, and the gas pressure is meticulously regulated at 4 
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Fig. 4. Model diagram of multiplicity devices established in MCNP.((a): XY cross-sectional view of MCNP model; (b): XZ cross-sectional 
view of MCNP model; (c): thermal neutron distribution map of XY cross section (d): thermal neutron distribution map of XZ cross section) 


atm. The summit plug comprises a stratum of graphite with a 
thickness of 15 cm. The lower region extends from the base 
of the sample chamber to the bottom of the device. The com- 
prehensive dimensions of the device encompass a height of 
82.7 cm and length and width dimensions of 74.2 cm. The 
model schematic is constructed using MCNP based on the 
aforementioned information and the device-design diagram 
in Fig. 1, as shown in Fig. 4 (a) and (b). 


Similarly, a *°*Cf source was positioned at the focal point 
within the confines of the sample chamber, and the dispersion 
of thermal neutrons on the axial and radial cross-sections at 
the central coordinates of the sample chamber was meticu- 
lously documented. Fig. 4 (c) and Fig. 4(d) illustrate the re- 
sults of the thermal neutron distribution at various positions 
within the device, as simulated using MCNP. 


Referring to Fig. 4 (c) and (d), because of the presence of 
moderators, fast neutrons rapidly lose energy and transform 
into thermal neutrons. Owing to the absorption of thermal 
neutrons by He gas in the neutron detectors, the thermal 
neutrons in the device achieve a sufficiently low level to be 
deemed inconsequential after they traverse the ultimate ring 
of the detectors. The incorporation of the Cd layer renders 
thermal neutrons within the sample chamber practically neg- 
ligible. Simultaneously, the Cd veneer on the outer surface 
of the initially decelerating structure absorbs residual thermal 
neutrons. In the outermost tier of the apparatus, thermal neu- 


trons are nearly nonexistent, thereby attesting to the extraor- 
dinary proficiency of the device in neutron shielding. 

The simulation conducted a preliminary calculation of the 
device detection efficiency, yielding a value of 42.10 % based 
on Eq. 15. According to Eq. 17, the detector die-away time 
of the device was computed as 54.66 us. These data serve as 
references for the subsequent calibration of the device. 


B. Passive neutron multiplicity device construction 


Upon determining the optimal parameters of the device, the 
device was constructed and assembled. The overall structure 
of the device is shown in Fig. 5 (a), and Fig. 5 (b) shows the 
upper computer-analysis software that supports the multiplic- 
ity device. Fig. 5 (c) shows the internal electronic circuitry of 
the device, where each *He neutron detector in FH-NCM/S1 
is connected to a preamplifier, totaling 80 channels positioned 
directly above the *He neutron detector. Three circuit boards 
are employed to collect information from the 80 channels of 
the *He neutron detector, and the circuit board is illustrated in 
Fig 5 (e). Each channel is powered independently to simplify 
maintenance and replacement. 

The signal-acquisition and transmission module comprises 
several key components: a He neutron detector unit, detector 
power supply, signal-adapter board (Fig. 5 (e)), and signal- 
acquisition board (Fig. 5 (d)). The *He neutron detector unit 
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Fig. 5. Overall and internal structural details of FH-NCM/S1 and the upper computer software. ((a) Overall appearance of the device; 
(b) upper computer-software measurement interface; (c) device electronics-circuit part; (d) FPGA circuit board;(e) signal-adaptation circuit 
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requires a power-supply voltage of 12 Vdc, which is achieved 
using an independent high-voltage power-supply module. In 
addition, the unit incorporates signal amplification, shaping, 
and a comparator output circuit. The neutron-detection sig- 
nal undergoes shaping and is output as a TTL-level negative 
logic signal, ensuring swift equipment response, even over 
long transmission lines. 


The chassis incorporates three signal-adaptation circuit 
boards, as shown in Fig. 5 (e). Each signal-adaptation cir- 
cuit board can provide 12 Vdc @ 500 mA output power, sup- 
porting up to 30 *He neutron-detector power channels. They 
are linked to the FPGA through a 40-pin terminal to scan and 
collect the output signals. 


The signal-acquisition board receives signals from the 
three-channel adapter board, overseeing the input voltage of 
the adapter board and the 12 V power supply voltage of 
the nine-channel *He neutron detector group. The interface- 
conversion chip transforms the pulse signals from the 80- 
channel *He neutron detectors, which are then collected by 
the FPGA. The FPGA concurrently monitors and conducts 
analytical calculations for the counting status of the 80- 
channel *He neutron detectors. The processed data output 
results are transmitted to the host-computer software through 
a field bus or 10/100 Mbps network port. Following the pro- 
cessing and calculation of the data using the upper computer 
software, the results are displayed in real time on the software 
interface, as illustrated in Fig. 5 (b). 


IV. DISCUSSION ON DEVICE CALIBRATION AND 
MEASUREMENT RESULTS 


A. Calibration Results of FH-NCM/S1 


I. Detection Efficiency 


Suitable measurement parameters are configured for the 
actual tests. The neutron yield of the **Cf source utilized 
in the measurement is 2.78 10° st. The source, including 
the background and known neutron yield, is measured for 
300 s six times. The total count of the *He detection array 
is recorded for each occasion, and the overall detection effi- 
ciency of the device is calculated using Eq. 15 and Eq. 16. 


After constructing the device, appropriate measurement 
parameters are set to consecutively measure the background 
and *°*Cf source with a known neutron yield six times, each 
lasting 300 s. The counting rate of the *He detection array is 
recorded for each occasion. The measured No is 1197267.3 
and N, is 1.5. The average count rate after deducting the 
background is 1.20x10° s!. Using Eq. 15 and Eq. 16, the 
detection efficiency (€) is calculated to be 43.07 %, which is 
close to 42.10 %, which is obtained from the simulation. 


2. Detector Die-away Time 


In the actual test, each measurement lasts for 20 s, and the 
measurement is repeated 30 times. The measurement results 
of the average value of S and double-rate D under different 
gate widths are recorded. Eq. 17 is used to fit the data ob- 
tained by using a curve-fitting method. The curve of Eq. 17 
is plotted, as shown in Fig. 7 (a). The fitting result for the 
value is 55.79 us, which is close to 54.66 us, obtained from 
the simulation. 


3. Sample Quantitative Results and Relative Error 


Owing to the insufficient availability of 7°?Cf sources for 
the measurement, the feasibility of the device was directly 
verified by comparing the measurement results of nine dif- 
ferent standard plutonium oxide samples (with *4°Pu masses 
ranging from 0.1g to 54.5 g) using the FH-NCM/S1 device 
and PSMC combined with JSR-15. These cylindrical plu- 
tonium oxide samples have the same density and encapsula- 
tion method, and the abundances of the different isotopes are 
consistent, with differences only in the diameter and height. 
The shift-register analysis mode of the FH-NCM/S1 device 
is used for the measurement and compared with the measure- 
ment mode of JSR-15. Each sample is measured for 300 s ina 
single measurement, which is repeated three times. Standard 
plutonium samples are measured and quantitatively analyzed 
using FH-NCM/S1 and PSMC devices. 

The relationship between the double rates and equivalent 
mass of **°Pu is shown in Fig. 7 (b). When the measuring 
materials are the same, the fitting curve of the double rate of 
the new device exhibits a shape similar to that of the PSMC 
device. However, the two curves do not completely overlap. 
According to the definition of the double rate (refer to Eq. 5), 
this discrepancy arises because it is influenced not only by the 
measurement material, but also by the detection efficiency of 
the device (€) and double-gate fraction (fa). 

To further evaluate the effectiveness of the device, accord- 
ing to Eq. 5, after eliminating the impact of the detection ef- 
ficiency (€) and double-gate fraction (fa), the correct double 
rate (D,) obtained when two devices measure the same sam- 
ple can be approximated as only related to the sample itself 
(Eq. 21). Therefore, the corresponding measurement results 
for the two devices for the same sample should exhibit a pro- 
portional relationship (De psmc= De FH-NCM/S1)- 
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The pre-delay times for FH-NCM/S1 and PSMC are set to 
3 us, the gate width is set to 64 us [28], and the long delay 
time is set to 2ms. The double-gate fractions (fa) of FH- 
NCM/S1 are calculated using Eq. 18 and the two parameters 
required for Eq. 21 and the die-away time are listed in Table 
L. 


Table 1. Parameters for measuring FH-NCM/S1 and PSMC 


Parameters O -o 
FH/NCM-S1 PSMC 
Detection Efficiency 0.431 0.533 
Double-Gate Fraction 0.6467 0.6279 
Detector Die-away Time 55.79 pus 49.10 us 


Fig. 7 (c) shows that the linearity of the double-rate (D) 
measured by the two devices for the same sample is excep- 
tional, indicating that the new device is capable of detect- 
ing Pu. Three potential sources of error exist between the 
two devices. First, errors may exist in the detection effi- 
ciency (£) and double-gate fraction (fq) of the device used 
to correct De. Second, the influence of dead time during the 
sample-measurement process cannot be ignored when the de- 
vice measures a large sample. Lastly, measurement errors 
may occur, particularly when the sample has a low number of 
fission neutrons. 


Samples No.1, No.2, and No.3, with 7*°Pu effective masses 
of 0.54 g, 27.23 g, and 45.50 g, respectively, are measured us- 
ing a single measurement of 20 s and 30 cycles. After 30 
measurement cycles, the calculated masses (m,) of 7*°Pu in 
the samples are calculated. Based on Eq. 11, Eq. 12, and 
Eq. 13, the results of the quantitative and error calculations 
for the three samples with similar appearances and different 
qualities are presented in Table 2. As shown in Table 2, the 
relative error between the measurement and calculation of the 
device is less than 5%, which satisfies the measurement re- 
quirements of the sample. 


Table 2. Measurement results and relative error of samples 


Sample 1 Sample 2 Sample 3 


Measurement ID mele Al% mele Al% mele Al% 


1 44.860 1.40 27.920 2.58 0.531 1.65 
2 44.706 1.75 27.063 0.57 0.557 3.13 
3 45.864 0.80 26.371 3.11 0.547 1.24 
4 46.060 1.23 27.755 1.97 0.533 1.33 
5 45.752 0.55 26.828 1.44 0.562 4.11 


4. The Repeatability of the Measurement Results 


The effective mass of **°Pu, obtained by measuring the 
same plutonium sample, is shown in Fig. 7 (d) and Fig. 7 
(e) with a single measurement time of 20s (a small measure- 
ment cycle). The relative standard deviation (RSD) of the 
calculation results is determined based on Eq. 21. The RSD 
value of the small-mass sample (sample No.3) measured us- 
ing the device is 7.23 %, and the RSD value of the large-mass 
sample (sample No.1) is 9.87%. The repeatability of the 
20s measurement is good, thereby satisfying the production- 
measurement requirements. 
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Fig. 6. Actual measurement results of FH-NCM/S1 ((a) Die-away-time fitting curve; (b) measurement curves of two devices for the same 
sample; (c) comparison of measurement curves after removing the influence of parameters; (d) multiple measurement curves of sample 1; (e) 


multiple measurement curves of sample 3) . 


B. Comparison of different analysis methods 


To test the accuracy of the transmission information of 
the pulse-position timestamp mode of the device, both shift- 
register analysis and pulse-position timestamp modes are em- 
ployed to measure the same unknown sample. After measur- 
ing for 300 s in each mode, the measured S/D/T and normal- 
ized distributions of R+A and A are recorded. The partial 
fragments of the pulse-sequence time spectrum recovered us- 
ing this software are shown in Fig. 6. The minimum resolu- 
tion for At in this mode is 20 ns. 


The distribution results recovered from the pulse-position 
timestamp mode and shift-register analysis mode using the 
same measurement parameters (pre-delay time is 3 us, gate 
width is 64 us, and long-delay time is 2 ms) are shown in 
Table 3. For clarity, the data in Table 3 has been nor- 
malized. The S/D/T values obtained by the two modes 
are 3119.20, 230.00, and 19.99 (pulse-position timestamp 
mode) and 3117.80, 230.20, and 17.11 (shift-register analysis 
mode), in cps. 

The CORREL function is used to compute the correla- 
tion between the distribution results presented in Table 3, as 
expressed in Eq. 22, to evaluate the coherence of the mea- 
surement outcomes across dual modalities. Within Eq. 22, 
Comprehensive 18 utilized as a parameter to gauge the unifor- 
mity of the two measurement methodologies, featuring values 
within the -l1—1 range. A proximity to 1 signifies heightened 


Table 3. Distribution of normalized information of (R+A) and (A) 
collected by two methods 


Pulse position Shift Register 
n timestamp mode Analysis Mode 
Pl) Qim)  P2m) Q2) 


0 7.27E-01 7.80E-01 7.27E-01 7.80E-01 
1 2.21E-01 1.86E-01 2.22E-01 1.87E-01 
2 4.40E-02 2.90E-02 4.36E-02 2.83E-02 
3 6.35E-03 3.85E-03 6.09E-03 3.53E-03 
4 9.62E-04 7.69E-04 1.03E-03 7.70E-04 
5 1.92E-04 0.00E+00 0.00E+00 0.00E+00 


congruity in the measurement results between the two modal- 
ities. ‘a’ and ‘b’ serve as invariable coefficients symbolizing 
the weights attributed to the correlation of the R+A and A sta- 
tistical data distributions, both set to 0.5. Pn and Qn signify 
the normalized distribution probabilities of R+A and A under 
the two modalities, respectively, as listed in Table 3. 


Ccomprehensive = @ Correl(P1, P2) + b Correl(Q1, Q2) 
\_(P1 — P1)(P2 — P2) 
\(P1 — P1)2 X (P2 — P2)2 
5 (Q1 — Q1)(Q2 — Q2) 


2 (Q1 — Q1) X (Q2 — Q2)? 
(22) 
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Fig. 7. Actual measurement results of FH-NCM/S1 (pulse-position timestamp mode). 


According to Eq. 22, the calculated value for Ccomprehensive 
is 0.999, demonstrating a striking similarity between the dis- 
tributions obtained using the two methodologies. This con- 
firms the viability of multiplicity data measurements using 
the pulse-position timestamp method employed by the device. 
Moreover, adopting this approach to collect pulse signals fa- 
cilitates convenient storage for offline measurements, data re- 
search, and analysis, thereby significantly enhancing the effi- 
ciency of data storage and organization. 


V. CONCLUSION 


This studyinvestigates the intricacies of simulation opti- 
mization, device fabrication, parameter calibration, and the 
nuanced evaluation of a novel passive neutron multiplicity 
apparatus, FH-NCM/S1, which is based on the neutron multi- 
plicity principle. The refined apparatus comprises a quartet of 
concentric rings housing *He neutron detectors. Each of the 
four rings accommodates 19, 25, 18, and 18 detectors, which 
are meticulously arrayed symmetrically around the central 
cavity. By harnessing the computational process of the Monte 
Carlo program MCNP, the optimization process entails fine- 
tuning the effective length of the *He neutron detector and the 
dimensions of the graphite moderator. The outcome indicates 
that the optimal length for the sensor is 30 cm, and the detec- 
tion efficiency can exceed 40 % with a moderator thickness of 
15 cm. 

Upon completion of the device setup, calibration was per- 
formed using a °? Cf source to determine the detection effi- 
ciency and die-away time, resulting in 43.07 % and 55.79 us, 
respectively. Subsequently, the shift-register analysis mode 
of the device was employed to measure nine plutonium ox- 
ide standards. After comparing the obtained dual counting 
rate (D) with the results obtained using the PSMC and JSR- 
15 combination, a modified dual counting rate (De) was pro- 
posed, validating the rationality of the shift-register analy- 
sis mode data-processing and solution methods. The perfor- 
mance of the device was demonstrated by measuring three 


standard plutonium samples using FH-NCM/S1 with a single 
measurement time of 20s and fluctuations between consec- 
utive measurements below 10%. After 30 measurement cy- 
cles, the relative error for 7*°Pu mass measurement was less 
than 5%, which meets commercial-application conditions. 

Finally, by comparing the measurements of the same 
plutonium-containing sample obtained with FH-NCM/S1 op- 
erating in both the shift-register analysis mode and pulse- 
position timestamp mode, the consistency of the results be- 
tween the two operational modes was verified. This study 
provides valuable information for the subsequent develop- 
ment, design, and performance assessment of passive neu- 
tron multiplicity measurement devices and offers meaningful 
guidance for future endeavors in the field. 
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